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Oh 1 ABSTRACT 
O 

t We have measured the rest-frame A ~ 1500 A comoving specific luminosity density 

of star-forming galaxies at redshift 3.5 < z < 6.5 from deep images taken with the 
Hubble Space Telescope (HST) and the Advanced Camera for Surveys (ACS), obtained 
^ ■ as part of the Great Observatories Origins Deep Survey (GOODS). We used color 

selection criteria to construct samples of star-forming galaxies at redshifts z ~ 4, 5 
and 6, identified by the signature of the 912 A Lyman continuum discontinuity and 
Lya forest blanketing in their rest-frame UV colors (Lyman-break galaxies). The ACS 
samples cover ~ 0.09 square degree, and are also relatively deep, reaching between 0.2 
and 0.5 L|, depending on the redshift, where L3 is the characteristic UV luminosity 
of Lyman-break galaxies at z ~ 3. The specific luminosity density of Lyman-break 
galaxies appears to be nearly constant with redshift from ~ 3 to z ~ 6, although the 
measure at z ~ 6 remains relatively uncertain, because it depends on the accurate 
estimate of the faint counts of the z ~ 6 sample. If Lyman-break galaxies are fair 
tracers of the cosmic star formation activity, our results suggest that at z ~ 6 the 



1 Based on observations obtained with the NASA/ESA Hubble Space Telescope obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under 
NASA contract NAS 5-26555. 
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universe was already producing stars as vigorously as it did near its maximum several 
Gyr later, at 1 z £ 3. Thus, the onset of large-scale star formation in the universe 
is to be sought at around z ~ 6 or higher, namely at less than ~ 7% of the current 
cosmic age. 

Subject headings: cosmology: observations — galaxies: formation — galaxies: evolution 
— galaxies: distances and redshifts 

1. INTRODUCTION 

The amount of star formation that took place early in the cosmic history, e.g. at z > 2, and 
its evolution with redshift are key pieces of information to constrain theories of galaxy evolution. 
Evidence of the possible signature of the trailing edge of the cosmic reionization epoch at redshifts 
z ^ 6 (e.g. Fan et al. 2002) has also recently renewed interest in measuring the star formation 
activity at these epochs, since starburst galaxies can contribute a significant fraction of the 
ionizing photons (Madau et al. 1999; Haiman et al. 2001; Steidel et al. 2001). Initial estimates 
based on the evolution of the cosmic neutral hydrogen mass density estimated from damped 
Lya absorption systems (DLA, Fall et al. 1996) suggested that star formation activity peaked 
in the past at around z ~ 1 and then decreased at higher redshifts. The first direct measures 
based on the observed UV radiation of relatively unobscured star-forming galaxies in the local 
and intermediate-redshift universe and of Lyman-break galaxies (LBGs) at z > 2 (Madau et 
al. 1996) seemed to confirm this general scenario. However, subsequent measures based on the 
systematic spectroscopic identification of hundreds of Lyman-break galaxies at z ~ 3 and z ~ 4 
showed that the star formation density traced by these sources, after reaching a maximum value 
somewhere between z ~ 1 and ~ 2, remains nearly constant up to the highest redshift observed 
with some confidence, i.e. z ~ 4 (Steidel et al. 1999, S99 hereafter). More recent measure at z ~ 5 
from deep, wide-area samples obtained with the Subaru telescope (Iwata et al. 2003) and with 
HST (Bouwens et al. 2003) show very mild evolution from z ~ 3, although possible evidence to 
the contrary has also been reported (Stanway et al. 2003a,b). Measures based on photometric 
redshifts also seem to be in general agreement with a constant or mildly decreasing activity of star 
formation as traced by starburst galaxies up to z ~ 6 (Thompson et al. 2001; Thompson 2003; 
Fontana et al. 2003; Kashikawa et al. 2003), or even suggest that the star-formation activity 
might actually increase with redshift (Lanzetta et al. 2002). 

The current measures of UV specific luminosity density of LBGs come from the Hubble Deep 
Field survey (HDF, Williams et al. 1996, 2000) and other deep, but relatively small survey made 
with HST (Bouwens et al. 2003) or from ground-based surveys. The HST surveys are limited by 
its very small area, and thus cosmic volume, which makes it prone to systematic bias induced 
by cosmic variance (for example, S99 argue that the original result by Madau et al. (1996) was, 
at least in part, the effect of cosmic fluctuations). The ground-based surveys cover a much 
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larger area, but at high redshifts the brightness of the night sky limits the sensitivity to only the 
bright end of the luminosity function. Unfortunately, LBGs are rather faint, with a characteristic 
magnitude z| 50 = 26.0 at z = 6 (based on Adelberger & Steidel 2000's UV luminosity function at 
z = 3) if no luminosity evolution occurs between z ~ 3 and z ~ 6. Thus, one needs to reach rather 
faint flux limits, even for 8-m telescopes, to account for most of the light from LBGs once the 
redshift exceeds z ~ 4. 

The GOODS/ ACS observations offer a good compromise, covering an area ~ 33 times larger 
than the combined HDFs, with much deeper sensitivity out to wavelengths of nearly 1 fim than 
that of most ground-based samples. This makes it possible to detect galaxies with luminosity 
fainter than L3 up to z ~ 6.5. In this letter we present measures of the comoving specific 
luminosity density over the range 3.5 z 6.5 based on the stack of three epochs of observations 
in both GOODS fields, the HDF-N and CDF-S. A companion letter by Dickinson et al. (2003, 
D03 hereafter) in this volume analyzes in greater detail the selection of galaxies at z ~ 6, including 
a discussion of the first spectroscopic identifications. Throughout, magnitudes are in the AB scale 
(Oke 1977), and the world model, when needed, is that of a flat universe with density parameters 
Q m = 0.3, = 0.7 and Hubble constant H = 70 km s _1 Mpc -1 . 

2. DATA AND SAMPLE SELECTION 

The data used in this letter are comprised of the first three epochs of observations, out of a 
total of five, of the GOODS/ACS survey (Giavalisco et al. 2003). They consist of eight mosaics, 
four in each of the two survey fields, and one in each of the -B435, Vqoq, 2775 and Z850 bands. Each 
mosaic covers an area approximately 10'xl6' in size with a scale of 0.05 arcsec/pixel. Table 1 of 
Giavalisco et al. (2003) lists the relevant parameters of the data, including the sensitivity. Source 
catalogs have been extracted using the SExtractor software, and the procedures are also described 
in Giavalisco et al. (2003). The catalogs used to select high-redshift galaxies have been created 
performing the source detection in the zs^q band and then using the isophotes defined during this 
process as fixed apertures for photometry in the other bands. We derive colors from isophotal 
magnitudes, but use the AUTO aperture-corrected magnitudes (in the SExtractor parlance) to 
measure the Z850~band apparent magnitudes. 

Samples of star-forming galaxies at high redshifts have been extracted from the matched 
catalogs using the Lyman-break technique (Madau et al. 1996; Dickinson 1998; S99, Steidel et al. 
2003; see also Giavalisco 2002). Specifically, we have defined LBGs at z ~ 4 (5435-band dropouts) 
using the colors equations: 

(B450-V6O6) > 1.2 + 1.4 x (V m6 -z 850 ) A (S 450 - U 606 ) > 1-2 A (U 606 -^85o) < 1-2 (1.1), 
and LBGs at z ~ 4 (V606~band dropouts) by: 



[(^606 - «775) > 1.5 + 0.9 X (i 775 - Z 8 5o)} V 
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V [(^606 - i 7 75) > 2.0] A (V 606 - i 7 75) > 1.2 A (i 775 - z 85 o) < 1.3 (1.2), 

where V and A are the logical OR and AND operators, respectively. For bands which are entirely 
shortward of the Lyman discontinuity (e.g. the -B435 for the Vejo6 _ an d ^75-band dropouts and the 
V(306 band for the i775-band dropouts) we also request a non detection with S/N < 2. 

These selection criteria are largely based on our previous experience with Lyman-break 
galaxies (e.g. Dickinson 1998; S99; Steidel et al. 2003). We have visually fine-tuned them based 
on the observed colors of stars and galaxies in the ACS images, as well as on the ACS synthetic 
photometry from galaxy spectral templates to reject most interlopers from lower redshifts, while 
efficiently detecting typical UV-bright, star-forming galaxies at redshifts of interest. 

Without a third band in the near-infrared, it is not possible to use a 2-color selection to 
define samples of LBGs at z ~ 6 in analogy to what done above. Instead, we have used the 
single-color threshold 

(i-z)>1.3, (1.3) 

along with a non-detection requirement S/N < 2 in the -B435 and V606 bands to define our sample 
of z ~ 6 ^775-band dropouts. D03 discuss the selection and robustness of these dropouts in detail, 
including the use of available near-IR photometry as "a posteriori test" of the selection criteria. 
A number of galaxies meeting these criteria have already been confirmed spectroscopically to be 
at z ~ 6 (Bunker et al. 2003; D03; Stanway et al. 2003b). Later on we discuss the implications of 
this different selection of LBGs on the measure of the specific luminosity density of the Z775~band 
dropouts. 

We have included in the samples only galaxies with S/N > 5 in the zs^o band, and we have 
visually inspected each of them and removed sources that were deemed as not real, either artifacts 
or spurious detections (estimated using counts of negative sources detected in the same data set). 
These amount to a negligible number for the -B435 and V6O6 samples, and ~ 12% for the ^775 
one. We have also eliminated from all samples sources with stellar morphology down to apparent 
magnitude ^850~ 26, i.e. where such a morphological classification is reliable. These amount to the 
3.1%, 8.3% and 4.6% for the -B435, Vqoq and ^775 samples, respectively. While the procedure bias 
our samples against LBGs that are unresolved by the ACS, it prevents contamination by galactic 
stars. In practice, it results in negligible changes to the measured specific luminosity density. 
Furthermore, while the vast majority of galaxies at z <C 6 have 2775-2775 < 1.3, photometric errors 
scatter some of these galaxies into our selection window. Thus, we have also purged the z'775-band 
sample from photometric-scatter interlopers using a statistical procedure. We have used a bright 
sub-sample of field galaxies, for which the photometry is very accurate, as a template of the color 
distribution of the sample as a whole and estimate the contamination using photometric errors as 
function of magnitude derived from the Monte Carlo simulations in each band (see D03). 

Down to ^85o< 26.5, roughly the 50%-completeness flux limit for unresolved sources (see 
Figure 4 of Giavalisco et al. 2003), the culled samples include 1115, 275 and 122 i?435-band, 
Vf306 _ band and z^-band dropouts, respectively. With a survey area of 316 arcmin 2 , this 
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corresponds to surface density S = 3.50 ± 0.10, 0.87 ± 0.05 and 0.39 ± 0.03 galaxies per arcmin 2 , 
respectively, where the error bars have been derived assuming Poisson fluctuations. 



3. THE UV COMOVING SPECIFIC LUMINOSITY DENSITY 

Because no systematic spectroscopic observations of the GOODS are available yet, we have 
used Monte Carlo simulations to estimate the redshift distribution function of the three samples. 
The technique consists of generating artificial LBGs distributed over a large redshift range (we used 
2.5 < z < 8) with assumed distribution functions of UV luminosity (we used a flat distribution, 
see later), SED, morphology and size. We have adjusted the input distribution functions of SED 
and size by requiring that the distribution functions recovered from the simulations using identical 
procedures as for the real Lyman-break galaxies match the observation at z ~ 4. In this way, both 
simulations and observations are subject to similar incompleteness, photometric errors in flux and 
color, blending, and other measurement errors. We have obtained the SED using the unobscured 
synthetic spectrum of a continuously star-forming galaxy with age 10 8 yr, Salpeter IMF and solar 
metallicity (Bruzual & Chariot 2000), and reddened it with the starburst extinction law (Calzetti 
2000) and E(B — V) randomly extracted from a gaussian distribution with He(B-V) = 0.15 and 
&E(B-V) = 0.15. For the cosmic opacity, we have used the Madau (1995) prescription extrapolated 
to z = 8. We have used an equal number of r 1 / 4 and exponential profiles with random orientation, 
and size extracted from a log-normal distribution function, as described by Ferguson et al. 2003, 
this volume, Figure 2). 

For a given dropout sample, the main output from the simulations is the probability function 
p(M, z, m) that a LBG with absolute magnitude M at redshift z is observed as having apparent 
magnitude m. A commonly used technique to derive the specific luminosity density of LBGs in 
this case is that of the "effective volume" (S99), where the spatial volume occupied by galaxies 
with apparent magnitude m in the sample is 

V eff (m) = J J p(M,z,m)dM^-^-dz. (2) 
The comoving specific luminosity density contributed by such galaxies is then estimated as 

d£(m) = n ^ L ^ ) dm, (3) 
V eff (m) 

where z is the average redshift of the simulated galaxies that have been selected into the sample, 
n(m) are the number counts of real Lyman-break galaxies observed with magnitude m, and 
L(m, z) is the specific luminosity of one such galaxy if placed at redshift z. The specific luminosity 
density is then L = J dC(m). The method provides a relatively accurate estimate in the case of 
the color-color selected S435~band and V606 _ band dropout samples, regardless of the assumptions 
about the UV luminosity distribution function used in the simulations, because these galaxies are 
selected based only on their observed SED, i.e. redshift, with little dependence of their intrinsic 
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luminosity. In other words, these LBGs have a relatively tight correlation of m and M up to the 
the detection limit. 

For the one-color selected 2775 -band dropouts the V e ff method underestimates the value of 
C, because in this case there is no tight correlation between absolute and apparent magnitude 
in the sample. Galaxies with a given apparent magnitude m now have absolute magnitude M 
distributed in a much larger interval, because whether or not they enter the sample depends on 
both their redshift and on their absolute magnitude. Specifically, the lower redshift bound is set 
by the color threshold, while the upper bound depends strongly on galaxy luminosity, as IGM 
opacity suppresses the Z85o~band flux at higher redshifts. Hence, the effective volume V e ff(m) 
is overestimated (unless in the simulation one uses the same intrinsic luminosity function fo the 
real galaxies), because in the simulations galaxies with widely different absolute magnitude are 
equally represented. In this case we have measured £ with a different procedure. We have done a 
X 2 minimization to find the intrinsic luminosity function cft(M) of the simulated galaxies such that 
once inserted into the real images and retrieved their number count 



is best-fit to the number count of the real galaxies. We use a Schechter function with slope fixed 
to the value found at z ~ 3 (Adelberger & Steidel 2000, AS00 hereafter), namely a = 1.6, and 
derive the parameter M* and (ft* from the fit. In practice, since our data only reach down to 
L ~ 0.2Lg at most, the assumption of a fixed value of a is a relatively minor one (we verified that 
values in the range 1.4 < a < 1.8 change our results by at most ~ 12%). The specific luminosity 
density is then simply C = J Lcft(M)dM, and for each redshift value this is computed using the 
corresponding sample's best-fit luminosity function (a detailed description of the fits and of the 
analysis of the result is the subject of a forthcoming paper). 

We found the average redshift and standard deviation of the redshift distribution to be 
zb = 3.78 and Sb = 0.34 for the 5435-dropout sample, zy = 4.92 and Sy = 0.33 for the Vejo6 
sample, and Zi = 5.74 and Si = 0.36 for the 1775 one. For each sample we have estimated C at 1500 
A (£1500) using K-corrections derived from a template LBG with the same UV color as the average 
one in the 5435-band dropout sample. We have derived £1500 integrating the luminosity function 
down to the "observed" limit, i.e. the absolute magnitude that for each sample corresponds to an 
apparent magnitude of Zg50=26.5, and we have also estimated the (relatively small) corrections so 
that all the measures reach the same absolute magnitude, which we have chosen to be that of the 
B-band dropout sample (0.2 L3). Finally, when using the V e ff method, we have also estimated 
the magnitude-dependent corrections for light losses due to finite aperture photometry that we 
have derived from the Monte Carlo simulations. 

The top panel of Figure 1 shows £1500 for the GOODS three samples derived using the two 
methods above, both the "observed" values and those corrected to the same absolute magnitude, 
compared to the measures at z ~ 3 and ~ 4 by S99. It can be seen that the two methods provide 
similar results for the S435 and V<306 dropout samples, but differ for the i775~band dropouts. 




(4) 
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Note that the luminosity function correction is relatively small for the V606 _ band dropouts, but 
increases for the i775~band ones, as the Z85o=26.5 limit of the sample corresponds to a brighter 
intrinsic luminosity at higher redshift. The bottom panel of the figure shows the star-formation 
density as a function of redshift. For the GOODS points we have used the conversion factor by 
Madau et al. (1998), namely SFR= 1.4 x 1CT 28 C in units of M yr" 1 Mpc -3 ; the other points 
are from S99, after conversion to our world model and correcting for the smaller range of absolute 
luminosity, i.e. > 0.2 Lg instead of > 0.1 L\. Both the observed values and those corrected for 
dust obscuration as suggested by AS00 are plotted. 

4. DISCUSSION AND CONCLUSIONS 

The specific luminosity density of LBGs appears to depend rather weakly on redshift over 
the range 2.5 ^ z 6.5. Integrating down to L ~ 0.2L|, S99 report Ci^qq(z = 3) = 1.50 ± 0.10, 
while from the GOODS data we find A500 ( z = 4 ) = 163 ± °- 05 > £i50o(^ = 5 ) = 1M ± °- 08 > 
and Ci5oo(z = 6) = 1.15lg!i9 m units of xlO 26 erg s" 1 Hz -1 Mpc -3 (1-cr error bars). However, 
while the points at z < 5 are relatively robust, the constraint at z ~ 6 is still somewhat weak, 
because it critically depends on the measure of the number counts near the sensitivity limit of the 
survey (see D03). For example, if we restrict the sample to Z850 _ band photometry with S/N > 6.5 
(roughly Z85o< 26) then we found that the best-fit specific luminosity density drops down to 
£ 1500 (z = 6) = 0.42t°; 33 x 10 26 erg s -1 Hz" 1 Mpc" 3 (0.42l°j° 2-a). At face value this is still a 
relatively mild drop from the z ~ 3 value, a factor ss 3.5, smaller than the factor f» 7 proposed by 
Stanway et al. (2003), but the error is too large for the constraint to be meaningful. We note, 
though, that the z ~ 5 point, which is more robust than the z ~ 6 one and does not need as big 
an extrapolation down the luminosity function, is reasonably consistent with the mild evolution 
scenario, lending support to our best estimate at z ~ 6. 

Thus, the result at z ~ 6 depends somewhat critically on the apparently substantial amount 
of light contributed by the faint galaxies. To the best of our knowledge the galaxies with 
5 < S/N < 6.5 are real z ~ 6 candidates. Based on our color selection, the observed dispersion 
of colors of faint galaxies, our simulations, and the expected colors of galaxies of various spectral 
types, we believe that our statistical corrections to the contamination is adequate and is not 
introducing a major systematic error in our measures. The very good agreement with the measure 
by S99 at z ~ 4, which is supported by systematic redshift identification of the galaxies, adds 
credence to this assertion. The current spectroscopic identifications at z ~ 6 (Bunker et al. 2003; 
D03; Stanway et al. 2003b) agree well with our predicted redshift distribution and also support 
this conclusion. The samples are still too sparse, however, to attempt a measure of the efficiency 
of the selection criteria and, thus, the contamination. Clearly, the measure of £1500 (z = 6) clearly 
needs to be revisited with deeper data. 

In any case, our measure at z ~ 6 is in overall good agreement with other similar measures 
from EST and ground-based data (Bouwens et al. 2003; Lehnert & Bremer 2003). Stanway et 
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al. (2003a,b) report a factor ss 7 decrease of the specific luminosity density, although this is also 
very likely in agreement with our resultof because they limit their measure to bright galaxies, i.e. 
L > L3. Actually, a direct quantitative comparison is difficult because of the difference in the 
sample selection and the pronounced dependence of £1500 on the faint counts. These authors also 
use the GOODS data, but they base their source detection on single-epoch images (we use a stack 
of three epochs), therefore necessarily reaching a shallower flux level and larger incompleteness 2 . 
Down to ^85o< 25.6 they find a total of 14 galaxies over 350 arcmin 2 (0.040 arcmin -2 ), while 
to the same flux level and using their same criteria we find 35 galaxies over 316 arcmin 2 (0.11 
arcmin -2 ). Also, note that they use the V e ff method to derive C, which, as we have detailed 
above, underestimates the specific luminosity density for the i775-band dropouts. 

Our measures are in good quantitative agreement with another deep, large-area LBG survey 
at z ~ 5 (Iwata et al. 2003) and also with surveys based on photometric redshifts (Thompson et 
al. 2001; Thompson 2003; Fontana et al. 2003; Kashikawa et al. 2003). They differ from some of 
the photometric redshift estimates (from the HDF) by Lanzetta et al. (2002), who, after correcting 
for bias in the photometry due to the (1 + z) 4 cosmological dimming, conclude that the specific 
star-formation density either remains constant up to z ~ 6, and then increases afterward, or 
increases monotonically. While we do explicitly account in our measure corrections for light losses 
due to finite aperture photometry (estimated from the Monte Carlo simulations), we do not find 
as large corrections as those suggested by Lanzetta et al (2002). Very likely, much of the difference 
is that our corrections are derived by requiring that the galaxy size distribution used as input 
in our Monte Carlo simulations is such that the output distribution matches the observations at 
z ~ 4 (this roughly corresponds to assuming that the sizes evolve as ~ H(z)~ 1 ; see Ferguson et 
al., this volume), while Lanzetta et al.'s corrections are based on the light distribution of galaxies 
atl^z^l.5, i.e. at considerably lower redshifts. 

It is not known if the UV specific luminosity density of Lyman-break galaxies at z > 2.5 is a 
fair tracer of the cosmic star formation activity (e.g. AS00). Regardless, if the dust obscuration 
properties of these sources are similar to local starburst galaxies (Meurer et al. 1999; Calzetti 
2000), and do not significantly change over the redshift range 3 < z < 6.5, then our result can be 
rephrased by saying that the star-formation activity of LBGs decreases very mildly with increasing 
redshift, at z ~ 6 being ~ 25% lower than it was around its maximum at 1 £ z 3, several Gyr 
earlier. Since LBGs are certainly sites of formation of a considerable amount of the present-day 
cosmic stellar mass density (AS00; Giavalisco 2002) it appears that the onset of substantial cosmic 
star formation takes place earlier than z > 6, namely at less than ~ 7% of the cosmic age. Finally, 
Figure 1 also shows that there is an overall good agreement between the dust-corrected data 
points and the semi-analytical models by Somerville et al. (2001, blue curve), and between the 
uncorrected data points and the predictions based on the observed evolution of the neutral HI 



2 They also use a slightly different color selection criterion, namely J775-Z850> 1-5. As D03 observe, this misses one 
of the galaxies spectroscopically confirmed at z = 5.8. 
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mass density as traced by DLA (which include the effect of dust obscuration) at least up z ~ 4, 
where these systems are known with goods statistics (Pei et al. 1999, red curves). 

What is the number of hydrogen-ionizing photons expected from our 2775 — band dropouts 
at z ~ 6? Since the rest-frame UV continuum at 1500 A is dominated by the same short-lived, 
massive stars that are responsible for the emission of photons shortward of the Lyman edge, the 
needed conversion factor, about one ionizing photon above 1 ryd for every 5 photons at 1500 A, is 
fairly insensitive to the assumed initial mass function and is independent of the galaxy history for 
ages 3> 10 7 yr (Madau, Haardt, & Rees 1999; Haiman, Abel, & Madau 2001). Here we normalize 
the number of ionizing photons to the observed 1500 A flux (rest-frame), i.e. we bypass the need 
for any correction due to dust extinction. A comoving luminosity density of £1500 = 10 26 er g s_1 
Hz -1 Mpc~ 3 implies then a comoving production rate of H -ionizing photons ^ion ~ 

3 x 10 51 s- 1 

Mpc~ 3 , or about 5/ csc photons per hydrogen atom per 5 x 10 8 yr escaping into the intergalactic 
medium. Here / e sc is the average escape fraction of ionizing radiation from the galaxy Hllayers 
relative to the escape fraction at 1500 A. Photoionization of intergalactic hydrogen requires more 
than one photon above 1 ryd per hydrogen atom, as extra photons are needed to keep the gas in 
overdense regions and filaments ionized against radiative recombinations. If / CS c is greater than a 
few tens of a percent, then our population of 2775 dropouts may contribute significantly to the UV 
metagalactic flux, and help to reionize the universe by redshift 6. 
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Redshift 




Fig. 1. — Top. The specific luminosity density at A = 1500 A of Lyman-break galaxies as a function 
of redshifts. The circles represent the 5435-band, V606 _ band and Z775-band dropout samples; the 
crosses are from Steidel et al. (1999). The hollow circles represent the V e ff measures, the filled 
circles represent the x 2 one - The blue symbols are the measures relative to the observed range of 
absolute luminosity; the red symbols include correction down to L = 0.2 x L| (see text). The error 
bars of the GOODS points are the 68% confidence interval of the fitting procedure. Bottom. The 
average star formation density of UV-bright star-forming galaxies as a function of redshift. The 
GOODS points have been obtained from the specific luminosity density using the conversion factor 
by Madau et al. (1998). The other points are from Steidel et al. (1999) after conversion to our world 
model. The red points are as observed, the blue points have been corrected for dust obscuration 
using the procedure proposed by AS00. The red curve is derived from the evolution of the HI mass 
density as traced by DLA absorbers (Pei et at. 1999; the blue curve is from semi-analytical models 
(Somerville et al. 2001). 



